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Abstract 

A litany of research has been published claiming strong solar influences on 
the Earth’s weather and climate. Much of this work includes documented 
errors and false-positives, yet is still frequently used to substantiate argu¬ 
ments of global warming denial. This manuscript reports on a recent study 
by Badruddin & Aslam (2014), hereafter BA14, which claimed a highly sig- 
nihcant {p = 1.4 x 10“^) relationship between extremes in the intensity of the 
Indian monsoon and the cosmic ray flux. They further speculated that the re¬ 
lationship they observed may apply across the entire tropical and sub-tropical 
belt, and be of global importance. However, their statistical analysis—and 
consequently their conclusions—were wrong. Specihcally, their error resulted 
from an assumption that their data’s underlying distribution was Gaussian. 
But, as demonstrated in this work, their data closely follow an ergodic chaotic 
distribution biased towards extreme values. From a probability density func¬ 
tion, calculated using a Monte Carlo sampling approach, I estimate the true 
signihcance of the BA14 samples to be p = 0.91. 
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1. Introduction 


Badruddin and Aslam (2014), hereafter BA14, recently reported a solar- 


terrestrial link between the cosmic ray (CR) flux and the Indian Monsoon, 
which they suggested may have implications of global importance and sup¬ 


port so-called ‘Cosmoclimatology’ (Svensmark, 2007). This work demon¬ 


strates the way in which their hndings were erroneous. 

BA14 based their claims on highly signihcant statistical relationships ob¬ 
tained from composite (epoch-superposed) samples. Specihcally, they exam¬ 
ined linear changes in monthly neutron monitor counts, analysed over m = 5 
month periods (during the months of May-September) from two samples, 
each comprised of n = 12 years of monthly resolution data: ie. composites 
from two matrices of n x m elements. The composites—which are vectors 
of the matrices averaged in the n-dimension—are respectively referred to as 
the ‘Drought’ and ‘Flood’ samples (which I shall also denote here as D and 
F), and represent the years of weakest and most intense monsoon precipita¬ 
tion respectively, recorded from 1964-2012. These data are shown in Figure 
[T| with the May-September periods of the composites emphasised: at hrst 
glance, it is true that these data show a linear change during the period 
highlighted, and also show anti-correlated between the D and F samples. 

Specihcally, BA14 evaluated the Pearson’s correlation coefficients (r-values) 
of the D and F samples, and used a standard two-tailed Student’s t-test 
(which assumes a Gaussian distribution) to test their probability (p) val¬ 
ues. They obtained values of r = —0.95 (p = 0.01) for D, and r = 0.99 
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Figure 1: Reproduction of the composite samples of BA14, showing the monthly-resolution 
pressure-adjusted neutron monitor count rate (units: counts min.“^ x 10^) from Oulu 
station (65.05° N, 25.47° E, 0.8 GV) occurring during 12 years of Indian monsoon ‘Drought’ 
(D) and ‘Flood’ (F) conditions. Composite means (in the matrix m-dimension) are plotted, 
with error ranges shown as ±1 standard error of the mean (SEM) value. The period of 
May-September, selected by BA14 for Pearson’s correlation analysis, has been emphasised 
in the plots. 

(p = 1 X 10“^) for F. Cumulatively, the p- value value was 1.4 x 10“^: i.e. 
such a result should occur by chance only 1/71942 times. BA14 interpreted 
these results to mean the lowest (and highest) precipitation volumes record¬ 
ing during the Indian monsoon period correspond to statistically signihcant 
decreases (and increases) in CR flux. 

From these apparently highly signihcant CR hux changes, BA14 con¬ 
cluded that a solar—monsoon link exists, and operates via a theoretical CR 
hux cloud connection. They speculated that this connection impacts the 
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monsoon in the following manner: increases in the CR flux enhance low 
cloud, rainfall, and surface evaporation, and also consequently decrease tem¬ 
perature (and vice versa). They further speculated that their hndings may 
be expanded to the whole tropical and sub-tropical belt, and as a result may 
impact temperatures at a global scale. However, the signihcance of the D 
and F samples—and consequently the conclusions—of BA14 are wrong. This 
error resulted from the assumption of a Gaussian data distribution, which is 
not true of their data, as I shall demonstrate. 

Moreover, of broader interest beyond the BA14 study is a recognition 
of a litany of fallacious solar—terrestrial studies: many of which have been 
re-examined in detail (e.g. by 


Pittock 

1978 

2009 

Farrar 

2000 

Kristjansson 


and Kristiansen 

2000 Damon and Taut, 

2004 

Sloan and Wolfendale, |2008 

Calogovic et ah 

2010 

Benestad and Schmidt 

2009 

Taken et al., 

2012bl 


Taken and Calogovic, 2013b). False-positives within this held are of partic¬ 


ular concern, as they contribute to a politically-motivated global warming 
denial movement. Providing material for groups intending to affect policy, 
such as the Heartland Institute’s Nongovernmental International Panel on 
Climate Change (NIPCC) or the Centre for Study of Carbon Dioxide and 


Global Change (Dunlap and McCright, 2010). Encouragingly though, a re¬ 
cent shift to open-access, and highly-repeatable workhows offers an oppor¬ 
tunity for rapid communal development (and cross-checking) across a broad 
range of helds, including solar-terrestrial studies: at minimum, such ap¬ 
proaches can more effectively facilitate the peer-review process, and enhance 
the quality and reliability of future publications. To illustrate this, this 


manuscript is supported by an accompanying iPython Notebook (Perez and 
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Granger, 2007), enabling users of the open-source software to easily check, 


repeat, and alter the analysis. This notebook (and all accompanying data) 


are openly available from hgshare (Laken, 2015). 


2. Analysis 


The CR flux oscillates between high and low values as solar activity pro¬ 
gresses from minimum to maximum during the ~ll-year solar Schwabe cycle. 
Consequently, over the 5-month timescales with which the BA14 study was 
concerned, the CR flux spends relatively little time at stable values. As a 
result, the population of r—values which can be derived from 5-month com¬ 
posites of these data are ergodic and biased towards extreme values. 

Using a Monte Carlo (MC) sampling approach, I have constructed 100,000 
composites of equal dimensions to the original D and F samples from the 
neutron monitor data, and obtained r-values over May-September periods: I 
note that the May-September restriction is not strictly required, as in reality 
the only requirement is that the MC-samples span an identical time-period 
(5-months) to the original samples. I refer to these data as Hq samples, as, 
by drawing these data randomly, they represent tests of the null hypothesis 
(for more details on this method applied to solar—terrestrial studies see 


Laken and Calogovic, 2013a). A probability density function (PDF) of these 
data are presented in Figure For comparison, a normalised Gaussian 
distribution—assumed by BA14—is also shown (dashed line). 

I have used two methods to model the PDF values: Firstly, a 4*^ order 
polynomial £t to the Hq samples (shown on Figure as the dotted line), 
of the function 0.03625a;^ - 0.0002797x3 - 0.0037^2 + 0.0007109x + 0.01964. 
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Figure 2: Probability density function (PDF) of r-values drawn from 100,000 null hypoth¬ 
esis (Hq) composites (from Monte Carlo samples of the Oulu neutron monitor data). For 
comparison, a normalised Gaussian distribution (with a mean and standard deviation of 
7.1 X 10“^ and 6.5 x 10“^) is plotted on the dashed line: BA14 wrongly assumed the 
data possessed this distribution, and consequently, this was the source of their error. A 
logistical map, which predicts chaotic distributions (given in Equation [^, is plotted on 
the solid line. A 4*^ order polynomial fit to the Hq population is plotted on the dotted 
line: this fit can be used to calculate the probability (p) of a given r-value. 
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And secondly, analytically using a Logistic map (as introduced by May et al. 


1976), which predicts the distribution of chaotically oscillating data (shown 


as 


the solid black line in Figure [^. The formula for this is given in Equation 
(Ruelle, 1989), where p is the probability, and u is the variable (in this case 
r-values). 


p{u) = 


7r\/T— 

The distribution of r-values appear to follow the Logistic map to a high- 
degree, indicating that the solar-cycle is oscillating chaotically. Indeed, the 


( 1 ) 


chaotic nature of the solar cycle has been well described (e.g. Mundt et al. 


1991 Kremliovsky, 1994 Rozelot, 1995 [Charbonneau 2001 [Hanslmeier and 


Brajsa, 2010 Hanslmeier et ah, 2013). Disagreement between the Logistic 


map and the PDF occurs at the most extreme values, where r<—0.9 or r>0.9. 


p=l- (0.03625a;^ - 0.0002797a;3 - 0.0037a;2 + 0.0007109a; + 0.01964) (2) 

As the polynomial ht accurately follows the PDF, it can be readily used to 
estimate the p-value associated with a given r-value using Equation From 
this, I calculate that the D and F samples possess p-values of 0.954 and 
0.948 respectively, resulting in a cumulative p-value of 0.91, i.e. a chance 
of occurring under the null hypothesis of 1/1.1. This result has a p-value 
four orders of magnitude larger than that estimated by the Student’s t-test 
approach of BA14, and is virtually guaranteed by chance. Consequently, 
I conclude that the high r-values obtained in the BA14 composites do not 
support a relationship between extremes in Indian precipitation during the 
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monsoon and co-temporal changes in the CR flux, but instead they are simply 
among the most commonly obtained values based on this sampling approach. 


3. Discussion 

The Cosmics Leaving Outdoor Droplets (CLOUD) experiment at CERN 
has demonstrated that ion-mediated nucleation may lead to enhancements 
in aerosol formation of 2-10 times neutral values under specific laboratory 
conditions—low temperatures characteristic of the upper-troposphere, and 
with low concentrations of amines and organic molecules—however, this ef¬ 
fect is absent under conditions more closely representing the lower tropo¬ 


sphere (Kirkby et ah, 2011 Almeida et al., 2013). Despite this, even if we 


assume that a significant nucleation of new aerosol particles form with so¬ 
lar activity, climate model experiments (which include aerosol microphysics 
schemes) have found that this would still not result in a significant change 
in either concentrations of cloud condensation nuclei or cloud properties. 
This is because the majority of the newly formed particles are effectively 


scavenged by pre-existing larger aerosols (Pierce and Adams, 2009 Snow- 


Kropla et ah, 2011 Dunne et ah, 2012 Yu et al., 2012[). These conclusions 


are supported by satellite and ground-based observations (e.g. Erlykin et al 


2009 Kulmala et ah, 2010 Laken et ah, 2012ai Benestad, 2013 Krissansen- 


Totton and Davies, 2013). For these and additional reasons, the IPCC AR5 


concluded that the CR flux has played no significant role in recent global 


warming (Boucher et ah, 2013) 


The numerous pitfalls into which solar—terrestrial studies in particular 


may fall, were lucidly outlined nearly 40-years ago by Pittock (1978). Despite 






















































































this, many studies with improper statistical methods, black-box approaches, 
and ad-hoc hypotheses still frequently appear. This problem is prominent 
within the held of solar—terrestrial studies. Consequently, the literature is 


replete with cases of demonstrated false-positives (e.g. Friis-Christensen and 


Lassen, 1991 Marsh and Svensmark, 2000 Shaviv and Veizer, 2003 Scafetta 


and West, 2008 Svensmark et ah, 2009 Dragic et ah, 2011), many of which 


have been (and continue to be) used as the basis for claims behind global 


warming denial (e.g. such as in Idso and Singer, 2009 Idso et ah, 2013), 
immediately making cases such as the one described in this manuscript a 
serious prospect in need of address. 
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